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Abstract In the Ising model on the simple cubic lattice, 
we describe the inverse temperature /? and other quantities 
relevant for the computation of critical quantities in terms 
of a dimensionless squared mass M. The critical behaviors 
of those quantities are represented by the linear differential 
equations with constant coefficients which are related to crit¬ 
ical exponents. We estimate the critical temperature and ex¬ 
ponents via an expansion in the inverse powers of the mass 
under the use of ^-expansion. The critical inverse tempera¬ 
ture /3c is estimated first in unbiased manner and then critical 
exponents are also estimated in biased and unbiased self- 
contained way including w, the correction-to-scaling expo¬ 
nent, V, rj and 7 . 

Keywords Ising model • expansion in mass • delta- 
expansion • critical quantities 
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1 Introduction 

Since the invention of lattice field theories, the border be¬ 
tween condensed matter models and field theoretic models 
is lost and the two fields of physics are unified Ifni2l . Tradi¬ 
tionary, the both systems are described in terms of /3 which 
indicates the inverse temperature or the inverse bare cou¬ 
pling constant. In the field theory side, we notice that the 
lattice spacing a or the equivalent full-mass can play the 
role of the basic parameter describing the models through 
the dependence of the bare coupling constant on a. This 
feature naturally appears in the large N limit of field theo¬ 
retic models 0 . Also for finite-A^ case of non-linear sigma 
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models, it was shown that the scaling behavior of the in¬ 
verse bare coupling has been captured from an expansion 
method in inverse powers of the full mass JU. Then, it is a 
natural step to take the reverse point of view into the con¬ 
densed matter models. That is, from mathematical analogy 
of field theory models, the inverse temperature /3 may be 
effectively described in the correlation length or the corre¬ 
sponding dimension-less mass parameter. The behaviors of 
the inverse temperature and other thermodynamic quantities 
(the magnetic susceptibility etc) would be described by the 
mass. The critical point is universally specified by the mass¬ 
less point and the scaling behavior described in the mass 
may emerge as the asymptotic behavior in the vicinity of the 
point. 

In the present paper, we consider the 3-dimensional (3D) 
Ising model on the simple cubic lattice and try to estimate 
critical quantities, the critical temperature and critical expo¬ 
nents by the use of sole information in the high temperature 
expansion which is equivalent with the expansion in the in¬ 
verse powers of the mass appropriately chosen. To the best 
of our knowledge, most of the unbiased and self-contained 
computation has been carried out in field theoretic models 
and computer simulations. For instance in Ising universal¬ 
ity class at three dimension, there exist only a few litera¬ 
tures 1511^ where pure series expansions (high temperature 
expansions without universality hypothesis) of Ising models 
themselves were employed for the full exponents computa¬ 
tion in a self-contained manner. As the unbiased analysis for 
the first time, Guttmann studied the 19th order high temper¬ 
ature series and obtained the estimates of /3c, 7 and v on 
three 3D lattices Q by the use of differential approximants 
method Q. After a decade, using the same method, Butera 
and Comi investigated 21st order high temperature series for 
the N vector model and estimated those three quantities IQ . 
In general, accurate estimation of critical quantities becomes 
possible under the bias of the inverse critical temperature /3c 
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and/or the correction-to-scaling exponent 9 obtained in other 
methods. It is desirable, however, to equip necessary tools in 
one computation framework and we like to propose a self- 
contained approach to critical quantities in the expansion in 
a mass parameter. 

We adopt as the basic argument the second moment mass 
M defined through the susceptibility x = Thj < > 

and the second moment p, = p < soSj > where Sj de¬ 
notes the Ising spin at the site j = {ji, j 2 j ja)- In more detail, 
consider the expansion of the two point function < 

sqSj > at low momentum p. Using the expansion = 

1 — ipj — \{pjy -I- • ■ we obtain < s^Sj >= 

X - = (2£'xVM)/(2£'x//i + + 0{p‘^)) 


(D — 3). Thus the squared mass M is given by 
2Dx^^ 

p p ■ 


M = 


( 1 ) 


In the high temperature expansion, x and p can be expanded 
in powers of /3 and the inversion of M(/3) gives the series 
expansion of /3 in 1/M. The key ingredient is the function /3 
expressed by M. The analogous relation in field theories is 
the inverse bare coupling expressed by aA where A means a 
finite mass scale. The dimensionless mass M plays the role 
of the square of aA. 

From the form of the asymptotic expansion near the crit¬ 
ical point, we consider that the critical behaviors of /3(M) 
and also other quantities are endowed with linear differen¬ 
tial equations (LDE) with constant coefficients |0 . The LDE 
describing the scaling law in terms of M is not satisfied 
by (3{M) at large M. On the other hand our information 
available is the 1/M expansion of P{M) with a certain con¬ 
vergence radius. Thus, it is not allowed to use LDE to the 
1 /M series. However, the so-called (5-expansion is expected 
to change the status drastically Q. Let us denote a func¬ 
tion of interest be f{M) and its truncated 1/M expansion, 
fy{M). As well, f^{M) denotes the asymptotic expansion 
valid in the vicinity of the critical point M = 0. Then, it 
has been verified in some models that ^-expansion induces 
a transform of fy{M) to />(f), and />(t) is indeed effec¬ 
tive to investigate the critical behavior of f{M), which can 
be confirmed by the plots of />(f) and its derivatives. In 
other words, there is an overlapping between the effective 
regions of />(/) and /<(t), and the critical quantities may 
be extracted by the LDE expected to be satisfied by /> (f) 
at the matching region of /<(t) and />(/). Eor a detailed 
analysis on the 2D square Ising model, see ref. |0. After the 
confirmation of the scaling region realized in /> (t) when 
necessary, we estimate critical quantities /3c, v, uj (= 6/v), 
rj and 7 in a fully self-contained way. 

This paper is organized as follows: In section 2, to be 
self-contained, we briefly review the method of 6 expansion 
on the lattice. In section 3, we shall estimate /3c from an 
expansion in the mass M in an unbiased manner. Rough es¬ 
timation of u is also obtained. In section 4, we investigate 


in detail an improved estimate of critical exponents. First, 
we estimate u and 77 by making use of the result w ^ 0.84 
summarized in ref. Col. Then, unbiased and self-contained 
estimation of the exponents w, v, rj and 7 is attempted. The 
last section is devoted to concluding remarks. 


2 ^-expansion 


Let a given thermodynamic quantity as a function of M be 
f{M) and has the limit, limM-s-o f{M) = fc. Then, our 
task is to estimate fc via the 1/M expansion of / to the 


order N denoted by />Ar = X]n=o ^nM We change the 
variable from M to / by 

l-(5 


M = 


t 


(0 < (5 < 1 ), 


( 2 ) 


and view the function in space where the critical region 
is enlarged by the factor 1 — (5. Then, we expand the function 
/>Af((l — S)/t) in powers of 6. Typical term in />Ar 
becomes as 


i{n + 1 ) 


^ 2 . 


.(3) 


=e(l + nS + 

By truncating the (5-expansion at 5^ and putting (5 = 1, we 
obtain 

(L + n)! 


M 




nlLl 


-f 


Thus, the (5-expansion of />Ar to order L gives 


N 


f>N ^ 


n—0 


{L + nf. 
n\Ll ^ 


(4) 


(5) 


Empirically, it is found that the truncation order L gives best 
result when L depends on N such that 131191 

L = N -n. ( 6 ) 

Throughout this paper we employ the prescription (| 6 ]). Then 


M" 




with the binomial factor 


A^! 


r{N + i) 


i\{N-n)\ r{n + l)r[N -n+iy 


(7) 


(8) 


We notice that (^) = 1 and The ^-expansion 

keeps a constant term invariant. 

Thus, we define the (5-expansion of />Ar(M) by 


N 

f>N{M) a„ 

n—0 


r =■. DM[f>N{M)] =: />7v(/).(9) 


To summarize, the (5-expansion creates a new function f>N{t) 
associated with />Ar (M) by transforming the coefficient from 
a„ to an(^) which depends on the order N. The symbol 
D N denotes the order dependent transformation from fy^iM) 
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to f>N(t)- Since (^) —>• ^ (n:fixed, N —>■ oo), one might 
think that the series becomes ill-defined in the N ^ oo 
limit. However, it is shown in some simple models that there 
exists a region of the form ( 0 , to) = I where />Ar(t) devel¬ 
ops a plateau there and converges almost everywhere to fc 
as iV —>• oo This is conceivable, since putting 5=1 
means M = 0 for every t (see (|2])). We thus assume that, 
when limM->.o f{M) = fc exists. 


lim />iv(t) = fc ( 10 ) 

N —^OO 

over a certain region I. As explicitly studied in ref. ID, in 
one-dimensional Ising model, the above result is analytically 
confirmed. In the square Ising model, numerical analysis to 
large orders of 1 /M expansion has also served us the con¬ 
vincing result supporting (fTOl) . In the case of simple cubical 
lattice, we can confirm (fTOl i by the direct numerical exper¬ 
iment for relevant functions as shown in various plots (see, 
for example, the plots of Z?Ar[/3>Ar] = ;9>Ar in Fig. 1). The 
point is, in the f-space, the whole region of I becomes the 
scaling region as long as the order N is high enough. We 
emphasize that the effective region does not include the ori¬ 
gin t = 0 inside. 

In the presence of phase transition, we must deal with the 
case where the expansion of / around M = 0 is not regular. 
Then we consider how in such a case the (5-expansion affects 
the small M behavior of f{M), supposing that / behaves at 
small enough M as / ^ /c + /iM“i -f / 2 M “2 + • • • =: 
/< where 0 < q;i < a 2 < ■ • When M = — 5) 

is substituted into and ((1 — 6)/t)°‘ is expanded in 5, 
giving M“ = -f ( 2(5 -f • • •), a reasonable truncation 

protocol for the best matching with /<Ar(f) is not found on 
logical grounds yet. Here, we proceed by simply extending 
the factor (|7| by the extension of F functions, 

( 11 ) 


where 

(N\ r{N + i) 

V-(2y “ r{-a + i)r{N + a + i)' 


( 12 ) 


We thus define the ^-expansion of /< (M) is given by 


DM[f<{M)] = fc + h 



H-. 


(13) 


The generalized binomial factor (_^) vanishes when a = 
1, 2, 3, • ■ • and 


As the exponent a is large, the factor decreases faster. All 
these facts appeal the advantage of /jv(f) compared to 

In the present work, our task is summarized to estimate 
fc and ttfe (fc = 1, 2, • • •) from the known series (|3- In the 
process, we use derivatives of /a? (= />Ar or /< 7 v)- To this 
end we remark that 


,|d„[/„(M)1 = d[(1/M)^/„(M) 


(16) 


which states that Z)A/-operation and differentiation is corn- 
mutable. It is convenient to use the following abbreviate no¬ 
tation, 

fN{t) = fl,^\t), (fc = 0,l,2,3,---). (17) 

From (|3 and the property of the binomial factor, we obtain 

f>l = Nif>N - f>N-l). (18) 


This relates the derivative to the backward difference of the 
functions with respect to N. It means that the stationary 
point of /> 7 v sits on the curve of f>N-i, one lower or¬ 
der function. The second order derivative obeys similar con¬ 
straint, 

= N{N - l)(/>^ - 2/>jv_i + />jv_2) + fi^l.(19) 

Hence, when f^j^ = f^^j = 0 at f = f*, it follows that 
f>N = f>N-i = f>N -2 at t = t*. In general, ([T8]l, ([T9]l 
and their higher order followers state on the convergence 
issue of the sequence of estimates. That is, if a point sat¬ 
isfying approximately the stationarity to the higher orders, 
the estimates of fc by that stationary point ensures the ap¬ 
proximate coincidence of estimations at successive orders, 
fc{N), fc{N — 1 ), • • •, showing good convergence behav¬ 
ior. 


3 Estimation of (3c and u 

The Ising model on the simple cubic lattice is defined by the 
action 

sl = l, ( 20 ) 

<ij> 

where the spin sum is over all nearest neighbour pair < 
i, j > on the periodic lattice. Our purpose here is to estimate 
the inverse critical temperature (3c defined in our approach 
by 


Oa([M“]= 0, (a = l,2,3,---). (14) 

This means that the regular contributions in the scaling re¬ 
gion are essentially eliminated by the 6 expansion. For non¬ 
integer positive power of M, the factor decreases with the 
order and disappears in the N —>■ oo limit, 

/ N\ Af-“ 


ft = lim B(M). (21) 

Past studies up to 2002 year on the estimation of (3c 
and exponents are reviewed by Pelissetto and Vicari Col in 
comprehensive manner. Of course there are important new 
contributions after ref. Qol on the subject of accurate esti¬ 
mations of the critical quantities II 11 II 1 2l[T3l[T4l[T5l[l 611 1 7l[T8l 
[1911^ . The best estimations of (3c to our knowledge before 











4 


Hirofumi Yamada 


2002 are /3c = 0.22165459(10) ED, /3c = 0.2216595(15) where 
Ea, and the recent works added /3c = 0.2216545(1) and 1 

0.2216542(4) ED, Pc = 0.22165455(3) E2, Pc = 0.22165463^8^ M' 


(25) 


lITSl and Pc = 0.221652(2) ll 20 l . Newer estimations agree 
with each other to the first 6 decimal places and, in this work, 
we refer the range of Pc as the modest one. 

Pc = 0.22165 ~ 0.22166, ( 22 ) 

within which all recent estimates are included. The world 
average of v up to the review ifTOl is given by v = 0.6301(4). 
Newer estimates are = 0.63002(10) lITSl . v = 0.630(5) 
ESI, V = 0.6306(5) El, V = 0.63048(32) El, v = 
0.629(1) ESI, V = 0.62999(5) El and i/ = 0.6301(8) 
ESI- Now, though most accurate estimates come from m 
and ED, we refer 

v = 0.6301(4) 

summarized in m, since new results support ( 

3.1 Preliminary studies 

Our approach is based upon the high temperature expansion. 
The magnetic susceptibility and the second moment have 
been computed up to /3^® by Butera and Comi ll^ . From the 
definition of the second moment mass squared E) and the 
result reported in El, X> = l+6/3+30/3^+O(/3^) and /r = 
6/3+72/32+0(/33), we find 1 /M = p+6p^ + fp^+0{P'^). 
By the simple inversion of the series, we then obtain 


(23) 


/3> = a; — 6 a;^ H —~ 312a;'^ 

O 


12596 


- 21432a;® 


1330848 


a;^ - 1745344a;® 


148384348 


-X 


797787336 


.10 


-15857888272a;®^- 

11748802870160 
7 ^ 


17341288504 ,, 

-n-" 

2106367479672 


,13 


13 


.14 


263968267347944 


,15 


15 


-186504592354608a;^®- 

-21535692193295224a;^® 

4449606807205690200 

19 ^ 

12821205881021198992 


33924951987330804 


,17 


17 


19 


.20 


197756701920466780928 ni 

-7-" 

3442869826889278353376 . 

-n- 

80156432259652309452520 
23 ^ 

116948936021276297965072 


.23 


,24 


10946582972904015563857296 


x^® + 0(a;^®)(24) 


The result of the J-expansion to the order N is readily ob¬ 
tained by multiplying the binomial factor to the nth order 
coefficient giving 


N' 


N 

/3>Ar := 7lAr[/3>Ar] = ^ 6„( jf". (26) 

Now, we turn to argue the behavior of P near the critical 
point to prepare our estimation work. In the study of phase 
transition via series expansion technique, we rely upon heuris¬ 
tic assumption of the power law near transition point. In ref¬ 
erence to EH, we start the arguments by the consideration 
of the Wegner expansion ED of the correlation length ^ at 
temperatures near l//3c in terms of the reduced temperature 
r = 1 — PI Pc- Let the spectrum of powers of corrections be 
generated by the basic ones, 1 , 6 , Oi, 62 , ■ ■ ■ and their integer 
multiples. Then, 


^ ~ fT~ 


1 -I- |aT®(l -I- aiiT -I- ai 2 T^ -I- 


} 


} 


(27) 


-h- 


iR-, 


-|-a2r^®(l -I- 0214" + a22T^ H-) 4 

-|-|6t®^(1 -I- biiT + 6i2r^ -I-) 

(1 + b2iT + &22'r^ H-) 4 

• -I- |mt(1 -I- UiT -I- U2t'^ + ■ ■ olj 

where = const x t(1 -|- tit -|- r 2 T^ -!-•••) stands for 
the analytic back ground and 0 < 6 * < 0 i < 02 < • • ■ 
The inversion yields a series in ^// with the spectrum of 
exponents of the form scaled by i/, —l/i'—{m+nd+ni6i + 
^202 + ■ ■ p/n where m, n,nk{k = 1, 2 ,3, • • •) are all non¬ 
negative integers. Leading term is (£,/provided 0 < 
1 and the next term is (5//)“*-^“''®^^’^. Here, we pose the 
assumption that 0 i > 1 and then the third term proves to 
be = (C//)"^^'" (provided 0 > 1/2). Then, 

inversion of the above series to a few orders reads 

T = (C//)-'/" + ^(C//)-('+®)/" + ^(e//)-"/" + - • • -( 28 ) 

The correlation length is defined through the large sep¬ 
aration limit of the two point function < s^Sj > (|j| —S' 00 ) 
and therefore has relation with the mass M defined in E)- 
Actually, it is known that for the case when the site j is 
on one axis, M = 2(cosh,^“^ — 1) ll25l . Thus, ^ 
/\^i/ 2 (i _|_ 0{M)) near the critical point and ^ log M 
in the high temperature limit (M — cxd). Though the high 
temperature behavior may be ~ const x \og{const x 
M) for the general position of the site j (One can under¬ 
stand this by directly performing the high temperature ex¬ 
pansion of < SQSj > at the leading order Il26l l. the rela¬ 
tion near the critical point should be universally given by 


25 
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Fig. 1 (Color online) Plots of P^n{x) (black) and /3>iv(f) (red) at 
N = 24 (dashed), 25 (solid). It is shown that /3>iv indicates the crit¬ 
ical temperature /3c = 0.22165 • • • in the plateau region. In /3> Ar(f), 
the limit of effective region may be signaled by the point of least vari¬ 
ation and at that point the value of /3c is best estimated. 



temperature region is pushed away to the neighborhood of the origin 
t = 0 and the scaling region appears around t ~ 0.11 for all plotted 
derivatives. 


~ + 0{M)), since the result must become 

independent of the direction of the site j. Thus, near the 
massless limit, ~ + 0{M)) and the term 

of order appears in r. This term belongs to rather 

higher orders in our study and thus we arrive at 

r = + + -fiM^+■■■), (29) 

and 


/3 = /3c(l - r). (30) 

Having finished a survey on the series, let us start the 
argument by supposing the scaling behavior of inverse tem¬ 
perature in the power-like form. By the combination of ( |29] | 
and (l30l l. we write 


/3< = p^-Aix-P^ -A2 X-p^ - Asx-p^ - O{x-p*),0^) 
where 0 < pi < p 2 < Pa < • • • and 

p<= Pc- Mt-P^ - A^t-p^ - Ast-p^ - 0{t-P^), (32) 

where 


An — An 


N 

K-Pn 

We note that 


1 


Pi = 


P2 = 


1 + e 


2i^ ' 

Since (_^ ) —0 as TV 
certain region of t. 


P3 = 


l-fl 

2i^ 


(33) 


(34) 


oo, we may assume that over a 


lim /3<Ar = Pc- (35) 

N —¥oo 

The effects of (5-expansion are clearly shown in the plots 
of relevant functions. For an example, the plot of P> 25 {x) 
and /3>25(f) in Fig. 1 shows that /3>25(f) exhibits the trend 
of approaching to the limit Pc, while /3>25(a:) has no sign 
of the scaling, as it would. This numerical examination pro¬ 
vides a clear evidence that the 5-expansion gives effective 
function of t capturing the scaling behavior at non-large t. 


~(k) 

We also show in Fig. 2 the plots of /3>25 and P ^25 foJ" ^ = 
1^5. Also in the derivatives, it seems that the scaling be¬ 
haviors have emerged; After oscillations of a few or several 
times, the derivatives set in their scaling behaviors of mono- 
tonically approaching to zero. It however needs a detailed 
study to judge to what level or grade the observed behaviors 
can be classified. On this respects, logarithmic plots of the 
derivative P^^'> (t) (k = 1,2,3,4) help the assessment. The 
log-plots of derivatives using Py'^ {t) effective at small t in¬ 
cluding the square model are shown in Fig. 3. The shown 
behaviors should be compared with the following, 

log((-l)'=+i^l;'=^) = log(pjAi)-pi \ogt+Oit-P^+PA06) 

We note that the linear-like behavior in the logarithmic plot 
is necessary for a benchmark of the scaling. If the behavior 
of log((—is linear and the correction 0 (f“P 2 +pi ^ 
is negligible, then the leading terms are dominant and the 
scaling behavior is pure and high grade. In this case, the gra¬ 
dient yields a good approximation for —pi — —l/{2iy). If 
the behavior is linear-like but the correction is still affective, 
the scaling behavior is a mixture of the leading and correc¬ 
tion terms. 

The square model at 25th order exhibits clear and pure 
scalings for the first and second derivatives, as shown in Fig. 
3a (We have used an exact mass M = 2(cosh^”^ — 1) 
1251 where the correlation length is given by /3 as = 
— log tanh P — 2P |]27|). The points appropriate for the lin¬ 
earization of the first derivative are at f = 0.172647 • • • 
and 0.210147 • • •. The slope values at respective points are 
—0.515381 • • • and —0.508829 • • •, respectively. As for the 
second derivative, the linearization is suitable at f = 0.190590 
and the gradient takes the value —0.544483 • • •. As expected, 
the value of the gradient in the first derivative at larger t is 
most accurate to the exact value —1/{21/) = —0.5. The third 
order derivative shows rather mixed scaling and the fourth 
order, weak and obscure. Thus, the times of differentiation 
increases, the scaling behavior becomes less-qualified. 
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(b) 


Fig. 3 (Color online) Logarithmic Plots of (fc = 1 ~ 4) in (a) 
the square Ising model and (b) the cubic Ising model al N = 25. 


At 3-dimension, the scaling behavior is not so high graded 
compared to the 2-dimensional case. For instance, though 
the curve of the third order derivative must be in the lower 
region of the second order one, it is located at the upper re¬ 
gion of the second order one. Thus, the trend of the scalings 
are indeed implied by the linear-like behavior, but the level 
of the scaling is not so high except for the first order. The 
vertical intervals between the derivatives are narrower at 3D 
compared with the 2D case. This is due to the fact that the 
interval roughly given by — logpi is smaller for 3D model 
implying the closeness of pi = l/(2i/) to unity. 

We emphasize that, the J-expansion is still effective to 
uncover the scaling region within small t regime. In par¬ 
ticular for /3>iv, it is obtained as the integration of its first 
derivative and the scaling behavior is expected to be better. 
The problem is that the leading correction of order ) 

is not negligible in ,3>Ar yet and the estimate by the sta¬ 
tionary point in the plateau does not lead high accuracy for 
N < 25. Now then, our idea upon the estimation of j3c is to 
cancel dominant corrections to fie in (l32l) by incorporating 
the derivative and consider the combination. 


l+Pi 


/3<- 


dlogx- 

Then it holds approximately at large x (small M) that 
d 


i+pr 


dlogx- 


^ 13c- 


(37) 


(38) 


The discrepancy between both sides is of order 0{x~p^ ) and 
reduced to the second order correction. 

From ([Tol l again, for large N, would approximately 
obey dMt over a certain region that 


l+Pi^ 


d - 
d log t. 


P<N = l3c + 0{t-P^), 


(39) 


where := Z)jv[/3<(M)]. We notice that, as the ex¬ 

tension, one can readily build up iTth order LDE such that 

K 


n[ 




dlogt 


^<N = Pc + Oit-P^+^), 


(40) 


for K = 1, 2, 3, • • •. For the sake of notational simplicity, we 
denote 


:= 


K 

n[ 

k^l 




dlogt 


P< 


N- 


(41) 


Making use of ( |39] | or (l40l i. we carry out estimation of /3c 
and 1 / in the following subsections. 

Closing comment on the naturalness of the scheme based 
upon LDE; We note that OTI) can be written as /3< = Pc — 
y^ig-pi iogx _ logx-form urges us to con¬ 

sider that the series is an expansion in eigen functions of 
the linear operator Hfcbfc + d/dlogx]. When there exists a 
quantity in which some terms are multiplied by the logarith¬ 
mic factor of X giving log x x x~^’‘, it means that the eigen 
value A/c is doubly degenerated. The magnetic susceptibility 
and the specihe heat at 2D just correspond to such a case. 
Also in field theoretic models, ordinary quantum corrections 
come as logarithmic corrections. The approach explored in 
this work would be capable of handling these cases. 


3.2 Estimate of Pc and v with the 1st order LDE 

The starting point is (IWb which states that <^< = [1 -f 
{d/d log t)]P^N approaches as W —>• oo to Pc- Eirst, at hnite 
N, we employ the principle of minimum sensitivity (PMS) 
ll^ . The use of the principle is quite natural since the con¬ 
vergence region of the left-hand-side of (l39l l or (l40l) should 
appear as a plateau. Thus, the value of ~ = P^n + 

^ possible stationary point gives an estimation 
of Pc- To obtain numerical results, the value of pi = l/{2v) 
is necessary. Hence, to hx all unknown quantities Pc and pi, 
we extend PMS condition such that the stationarity of 
becomes maximal at a proper value of pi by further impos¬ 
ing the second derivative = P^]^ -f pj"^/3^^ to van¬ 
ish or be minimum in the magnitude (This extension was 
first proposed in ref. 1291 ). Now, the point is to remind that 
M exhibits the scaling behavior in some region. From Fig. 
1 , we convince ourselves that the upper limit of the region 
showing the scaling behavior of Py\^ {k = 0,1, 2) is about 
t ~ 0.15. Making use of this, we replace /3<Ar by Pn^ in 
the estimation task under extended PMS. 
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-(k) 

Explicitly the extended PMS condition in terms of ^ 
reads 


=^>]/+Pr'^>lt = 0, (42) 

(43) 


The symbol ~ in the second condition means the local min¬ 
imality in the magnitude. It is straightforward to obtain 
from (l42ll as the function of t which indicates the estima¬ 
tion point. To avoid possible confusion we denote the so¬ 
lution as p(t). Thus, p{t) = Then, substitut¬ 
ing it into (021, we obtain — 

P>n^>n)/^>n fro™ ^ lp(t) ~ 0 the point t = t* at 
which the stationarity is maximally realized is found. Next, 
p{t*) = 1 /pI gives the estimate of by i/ = l/(2pj). In 
this manner, we have estimation of /3c by 


/3c = '^> 


= + (1M)/3W (r 


(44) 


We carry out this procedure from 4th order to 25th order. 
The reliable and steady trend of our approach emerges from 
13th order. The trend is the characteristic scaling behavior 

—(o'] 

that the p(f)-substituted second derivative /3>jv -f p{t)j3/^jq 
should exhibit, which we describe below: The solution p{t) 
should behave at t in the critical region as 


pit) = 


Y.V=iiPk)^Aut-v^ 


1 

Pi 


P 2 A 2 

Pi^i 


iP2 


p^]t-iP2-pi) ^ _ 


(45) 


Then loses t term and has to 

behave as 

!^> ) = —(pi -P2)‘^A2t~P^ H-. (46) 

Pi 

Note that p(f) ^ 1/pi since (_pj/(_pj ~ r{l-p 2 )/r{l- 
Pi) X A^“P2-i-pi tends to zero (pi < P2) as — >■ cxd. In the 

same manner, the correction in (l46T l vanishes in the N ^ 00 
limit. At finite order, the correction ~ const x however 
remains. For instance, as implied in Fig. 4(a), the weak peak 
of + p/t/Pylf about t ~ 0.09 at = 13 is 

interpreted to show the transition from high temperature to 
the scaling region. The dotted vertical lines in Fig. 4(a) and 
(b) represent the divergences of and p(t) at N = 13. 

As shown in Fig. 4 (a), (b), we obtain two candidates 
of the /3c estimation. One comes from the zero of = 
Pylf3-pit)p)^lf and the other fromits local minimum (^ 0). 
The first candidate is outside (t ~ 0.07) and the other is at 
the beginning (t ^ 0.12) of the scaling region. The later 
one represented by the blue cross is also the point mark¬ 
ing the upper limit of the reliable region of the truncated 
pit) 'P>(f). Hence, we regard the candidate 
marked by the blue cross as important and optimal. Actual 
values computed are /3c = 0.218638 (at t* = 0.068789) 



.0.02 I ' ■ ' ' ' I. . I I I I . 

0.00 0.05 0.10 0.15 0.20 0.25 



(b) 



ft>Nl^>N)/l^>N’ (b) pit) and (c) ^y{p*,t) at N = 13. Vertical 
dashed lines in (a) and (b) indicate the singular points coming from the 
solution of = 0. They are at f = 0.01526 and 0.14768. Note 
that between the two points, p is positive definite. The range of optimal 
stationary points are thus limited within 0.01526 < t* < 0.14768. 
Two candidates of physical interest due to the extended PMS are indi¬ 
cated by the red triangle and blue cross. Among them, important one is 
the solution represented by blue cross, which we call proper solution, 
because it lies in the place where the critical region begins to appear. 
The lower plot shows two !P> (p* ,t) at N = 13. The red one is for 
p* = 0.863705 and the blue one is for p* = 0.945108. 


and /3c = 0.220024 (at t* = 0.116908). To summarize, we 
can thus select which solution is the best choice among ap¬ 
peared. We note that at the point where = 0, p{t) be¬ 
comes extremal. At general order N, this is found by the 
direct differentiation of p{t) by logf. Turning to the even 
high order case N = 24, from the plots shown in Fig. 5, we 
find that the point marked by the blue cross gives the best 
solution where =0. 
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Table 1 Estimation of /3c and pi = l/(2i') = 0.79365 (1/pi ~ 1.26) with one parameter. Only the proper values are shown. 


order 

20 

21 

22 

23 

24 

25 

/3c 

0.220933 

0.220944 

0.221035 

0.221040 

0.221114 

0.221117 

1/pi 

1.028607 

1.029841 

1.041745 

1.042425 

1.053098 

1.053408 

t* 

0.133502 

0.126393 

0.135082 

0.128043 

0.136505 

0.129586 



.0.021 I ' ' i.1. I 

0.00 0.05 0.10 0.15 0.20 0.25 

(a) 



(b) 


Fig. 5 (Color online) Plots of (a) = ((/3 >m)^ - 

(l>) p(i) at N = 24. Scaling regions have 
developed broader compared to the 13th order. This agrees with the 
assumption (H. 


We call the solution marked by the blue cross as the 
proper solution and take it as giving the legitimate estimate 
of f3c- In this manner, we can identify the proper sequence of 
estimated /3c. In Fig. 6, the sequence of /3c estimated at non¬ 
proper t is shown by red triangles and at proper t by blue 
crosses. It is clearly shown that the proper estimates serve 
us better values. 


0.224 

0.222 

0.220 

0.218 

0.216 

Fig. 6 (Color online) Estimated result of /3c from 4th to 25th orders. 
The channel of proper sequence opens at At = 8 for even orders and 
At = 13 for odd orders. 


A 
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20 


25 


rections are composed by first few or several terms in (l32l l 
and not only by Ait~P^. Hence, even when the cancelation 
is successful, it does not necessarily mean that 1 /pj contains 
accurate information of true 1/pi = 2t/, unless the order of 
1 /M expansion is extremely large. 

The sequence of good estimations starting with 8th or¬ 
der for even N and 13th for odd N has enough number of 
terms to do extrapolation to the infinite order by a fit assum¬ 
ing the simplest form,/3c (V) ^ /3c —Using the last 
three estimations of odd orders, 21st, 23rd and 25th, we ob¬ 
tain /3c = 0.221642 with b = 0.101171 and u = 1.634371. 
If we use the 20th, 22nd and 24th orders results, we obtain 
/3c = 0.221624 with b = 0.10283 and u = 1.669628. Both 
values come to close to the established value, /3c = 0.22165 
(see (I22I 1). The relative error is only 0.00035 ~ 0.01 %. 
Also for we have done the fitting and obtained with 
21st, 23rd and 25th order results, p/”^ ss 1.2828 and with 
20th, 22nd and 24th orders, p/”^ r; 1.2640. The estimation 
is much improved as having relative discrepancy to 1.2602 
about 0.3 ^ 1.8 %. We can say that the proper sequence 
tends to the present world standard. 


As the order increases, the optimized value of pi gradu¬ 
ally increases but is not sure whether the sequence converges 
to the established value, 1/pi — 21^= 1.2602. See Table 
1. On the other hand, the accuracy of estimated /3c is good 
for N = 25 as /3c = 0.221117 ■ • • and the relative etTor is 
about 0.24 %. The non-accurate results for pj"^ may be ex¬ 
plained as follows: In our method of using extended PMS, 
the value of pi is fixed to achieve the total corrections being 
suppressed. Here notice that the dominant part of total cor- 


3.3 Estimates of /3c and u with the 2nd order LDE 

In this subsection, we apply our method to the ansatz /3< = 
/3c + const X + const x t~P'^ which satisfies the fol- 


lowing LDE, 



t d 1 

r t d 1 

= 


dlogJ 
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“ ~(k) 

By using /3>Ar and its derivatives in the place of we 

must search three optimal values of pi, p 2 and t. For the 

purpose, we use extended PMS specified by the couple of 

equations, 

=P>l + P^% + <yP>N = 0, 

= P>N + pP^n + - 0 > ( 48 ) 

where 



. 0.03 .. . . . i .. 

0.00 0.05 0.10 0.15 0.20 0.25 


P=Pl^+P 2 \ {P1P2) ^ 


(49) 


It is interesting to see the results of estimation, when 
the 5th order derivative, that shows the behavior insufficient 
for the scaling, enters into the estimation program. As in the 
case of one-parameter ansatz, it proves convenient to express 
p and (T in terms of t. From the first two equations, we obtain 


P{i) 

(j{t) 


0(2) 0(3) _ 0(1) o(4) 
P>nP>N P>nPn 

A 

5(1) 5(3) _ ( 3 ( 2 ) (2 
P>Ny>N \P>n) 

A 


where 




(3) ^2 


(50) 

(51) 

(52) 


The two parameters behave at large t as 

p(t) ~ + const X f“P3+P2 -I- ... ^ 

P 1 P 2 

a{t) ~ — -h const X + .... (53) 

P 1 P 2 

The critical behavior of + + 

is then given by ^ const x f“P^ and decays faster than f“Pi . 
It should also be noted that the coefficient of fP^ denoted 
by ’’const” depends on the order N and tends to zero in the 
N —>■ 00 limit. Therefore the 3rd order derivative in which 
the solutions dsoli and dSB are substituted should exhibit 
rapid trend of vanishing. At N = 25, this is numerically 
confirmed as shown in Fig. 7. The vertical dotted lines indi¬ 
cate the singularities of the function. All of them come from 
zero points of A, the solutions of Z\ = 0. The proper station- 
ary point lies between the two singular points. Actually, ' 
stays close to the horizontal axis between the two marking 
points (red triangle and blue cross) and the expected scaling 
(though not so clear) is observed there. Thus, the proper so¬ 
lution comes from the point represented by the blue cross. 
Now the solution of t is obtained by searching zero or min¬ 
imum point of + p{t)P^^]^ + a{t)P^lf. Then, 

with obtained t*, p* = p(t*) and a* = CT(f*), the estimation 
of /3c is directly given by 


Pc = P> 


— P>N + P* 


(54) 


Fig. 7 (Color online) (p(i)> o'(t), t) at N = 25. In the plot, we 
emphasized two candidates both of which provide positive two param¬ 
eters p and (j. The blue cross at the local minimum point indicates the 
proper solution which is considered as lying in the scaling region. The 
red triangle indicates the other, non-proper solution. 



Fig. 8 (Color online) !P> (p*, cr*, f) for p* = 1.6896423, cr* = 
0.5988992 (blue curve) and for p* = 1.5529395, cr* = 0.4763825 
(red curve). The dashed line indicates /3c = 0.22165. 



Fig. 9 (Color online) Plots of estimated /3c with two-parameter ansatz. 
The dashed line indicates /3c = 0.22165. The blue cross indicates 
proper estimate which appears from 20th order. 


The result of estimation is shown in Table 2, Figs. 8 and 
9. The behavior of ' becomes steady at 12th order and 
single solution is obtained up to 19th order. The new chan¬ 
nel to the proper sequence opens from 20th for even orders 
and from 23rd for odd orders. The level of the scaling behav- 
ior observed in ’ is lower compared to the one-parameter 
ansatz. For example, the scaling level of 25th order in the 
two-parameter case is, to the eye, the same level with the 
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Table 2 Estimation of /3c, pi and p 2 with two parameters. 


order 

20 

21 

22 

23 

24 

25 

/3c 

0.221442 

0.22140 

0.221505 

0.221513 

0.221553 

0.221562 

1/Pi 

1.138882 

1.127291 

1.160289 

1.163069 

1.179823 

1.183677 

1/P2 

0.433345 

0.413636 

0.468705 

0.472926 

0.500057 

0.505965 

t* 

0.116906 

0.096589 

0.119071 

0.114190 

0.120643 

0.116991 


13th order in the one-parameter case. However, we notice 
that the accuracy is improved; The relative error of 25th or¬ 
der j3c is about 0.004 %. The estimation of = 2^ is not 
so good yet, though two-parameter estimation is improved 
compared to one-parameter case. At 25th order, the relative 
error is about 6 %. 

The extrapolation of the proper sequence to the N ^ oo 
limit is not adequate in the two-parameter case. This is be¬ 
cause the number of elements in the sequence is not enough. 


3.4 Estimates of /3c and v with the 3rd and higher order 
LDEs 


The three parameter ansatz takes the form/3< = Pc+constx 
+ const X + const x t~P^. With third order LDE, 


<F< 


n [i+p*' 



/3c, 


(55) 


we continue estimation of /3c and u. The estimation proce¬ 
dure follows those of 1st and 2nd order LDEs and we omit 
the details and present just the outline. The extended PMS 
condition reads 


= 0 , 

+ PP>N + 

+ PP>N + = 0 . 

+ P^>N + + '^P>N ~ 

where p, a and r are parametrised as 

1 1 1 

p = -1-1-, 

Pi P2 P3 

1 1 1 

(T —-1-1-, 

PlP2 P2P3 PSPI 

1 

T = -. 

P 1 P 2 P 3 


(56) 


(57) 


By the first three equations of (l56l l. the three parameters p, a 
and T are obtained as functions of t and then the substitution 
of them into the last equation of (l56l l gives the solution t = 
f*. Then p, cr and T are estimated by p* = p(t*),a* = a(t*) 
andr* = r(f*) and /3c is obtained by = /3>Af+ p*/3^^^ + 
cr* -I- T-* I 


To 22nd order, the fourth order derivative shows a 
complicated behavior. Only from 23rd order, our method be¬ 
gins to provide estimations characteristic to the 3-parameter 
ansatz. The behavior of is not developed yet compared 
to the higher order results in the 1- and 2-parameter cases. 

It is typically reflected to the narrowness of the region (p, 
a and t are positive there) where the scaling is expected to 
emerge there at higher orders. In addition, at the last even 
order N — 24, the behavior of resembles to the 12th 
or 14th orders in the two-parameter ansatz. See the plot (a) 
in Eig. 10. These features are understandable since three- 
parameter ansatz involves 7th order derivative (see (|56|) which 
shows the behavior far from the scaling to 25th order. 



.0.10 I - . I. J J. !, , I I . 
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(a) 



= 0 indicated by the red triangle, t* = 0.1129729, is the 
only one (though not the proper one yet) to be considered as represent¬ 
ing the legitimate estimation. The solution gives p* = 2.1339308, 
a* = 1.2321906 and r* = 0.1743783. (b): {p*, a*, t* , t) is 

plotted at A = 25. The red triangle gives estimate of /3c = 0.22165, 
which is indicated by the dashed line. 
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Table 3 shows the results at 23rd, 24th and 25th orders. 
These three estimations exceed j3c = 0.22166, though the 
last order estimation is most close to the established value. 
We notice that at 25th order, the estimation of is also 
most accurate to give v = 0.6372. This implies that when 
Pc is precisely obtained, estimated v is also accurate. The 
opening of the proper and accurate estimate for the 3rd order 
LDE case demands further computation of high temperature 
expansion, maybe up to 30th order or more. 

Under the simple estimation without using extrapolation 
to the infinite order, increasing the number of exponents in¬ 
corporated improves the estimation of critical quantities so 
far. However, the reliable estimation with confidence of the 
scaling behavior of relevant functions sets in at larger orders 
when the number of exponents in the ansatz are increased. 
This stems from the fact that the differentiation on PN{t) 
creates oscillation and delays the appearance of the scaling 
behavior. Let us briefly argue the points: The small t behav¬ 
ior of is given by P^\x) = (—-|- 

P2A2t~P^ H- ). From 0 < pi < p2 < ps < ■ ■ ■ and the re¬ 

sult in Table 3, we find p 2 > 1. Then for n>2,p^{n> 2) 
grows with k. This means that the differentiation enhances 
higher order corrections and the critical behavior of /3< is 
obscured. While at small t, the differentiation on creates 
n to the coefficient of f”. Hence, also in small-f expansion, 

~(k) 

the upper limit of effective region of ^ tends to shrink. 
We have actually found that, up to 25th order, m-parameter 
extension when m > 4 does not work (In the 4-exponents 
ansatz, p^^'> enters into the estimation task and the function 
does not show any sign of the scaling to 25th order). In our 
study up to 25th order, the 3rd order LDE with three expo¬ 
nents is at the limit of our method. 

It would be better to mention on the p2 estimation. As 
would be understood from Table 2 and Table 3, the estima¬ 
tion of p2 is not successful yet. If one uses standard values 
for 1 / = 0.6301 and 6 — 0.5, one has = 2v j (1-1-0) ~ 
0.84. Our estimated results in two- and three- parameters 
ansatze are still far from the value. In the following section, 
we turn to the improved estimation on the critical exponents 
including the correction to the scaling w = 9/v. 


Table 3 Estimation of Pc and pi and p 2 with three parameters. 


order 

23 

24 

25 

Pc 

0.222079 

0.221741 

0.221687 

1/pl 

2.010008 

1.336976 

1.274452 

1/P2 

0.920863 

0.715954 

0.648486 

1/P3 

0.299087 

0.242590 

0.210993 

t* 

0.111805 

0.111116 

0.112973 


4 Improved estimation of critical exponents 

We here attempt to improve estimations of critical expo¬ 
nents. The reference values are (l23t and the following sum¬ 
marized in ref. Col, 

w = 0.84(4), (58) 

p = 0.0364(5), (59) 

7 = 1.2372(5). (60) 

The number in the parenthesis indicates the uncertainty or 
error in the last digit. Recent works added further estima¬ 
tions; uj = 0.782(5), T] = 0.0318(3), 7 = 1.2411(6) from 
m, w = 0.832(6), 7 = 0.03627(10) from Ha, w = 
0.82(4), 7 = 0.034(5) from Qa, w = 0.80(1), rj = 0.0341(5) 
from na, a; = 0.8303(18), 7 = 0.03631(3) from HD and 
oj = 0.83(9) from ll20ll . Though some results are out of the 
range indicated in (fSST l to (I6OI 1. we use those values quoted 
in nni as the bench mark. 

4.1 Preliminary studies 

From the lesson of the previous section, we learned that the 
massless limit of a given function itself is approximated bet¬ 
ter than the exponents of the corrections included. Thus, we 
use functions derived from P{M) and xi^) which provide 
associated exponents as the leading terms in the M —> 0 
limit. 

To improve the estimation of the critical exponents, we 
will use the characteristic structure of (I30ll and LDE, 
leading self-consistent point of view for v. To prepare the 
estimation work, we study the expansion structure of rele¬ 
vant functions appropriate for the estimation of critical ex¬ 
ponents. 

For the estimation of p, the ratio /P^^'^ =: fp is 

convenient since i/ itself appears as the leading term, 

= fp< = + ■ ■ ■ • (61) 

The amplitudes Bi, B 2 , • • • are written by a, v, 6 and so on 
(see (I27I 1') but details are not relevant for our purpose. 

The magnetic susceptibility x is defined by 

X = ^ < SoSn > ■ (62) 

n:sites 

To obtain x in the critical region expressed in terms of x, it 
suffices to substitute t{x) in (|29]) into the standard expres¬ 
sion of x<, X<(t) ~ C't”'>'[1 -I- const{T^ -|- • • •) -f ■ • • -f 
const{T -1- • • ■)] -I- xr iXR denotes the analytic back ground) 
IMI . The result reads 

7 0 1 

X< = Cx^ (1 -I- const ■ x~^ + const ■ x~^ -!-■••). (63) 

Scaling relation due to Fisher m reads that 7/(2z/) = 1 — 
ry/2 and what we can directly estimate is p rather than 7. In 
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the estimation of ?7, it proves convenient to address (logx)^^^ = 
{d/d log x) log X ■= fx- behaves near the critical point, 

fx< = ^+ + ■■■. (64) 

All behaviors of the series (IhTT l and (l64l l can be written 
as 


f<= fc + fix + f2X ® H-, 

where 

w 1 

The function /< satisfies ATth order LDE, 

K 


:= 


n[ 

n=l 


1 + g. 


-1 


dloga;. 


/< = /c + 0{x 


-gjf+i' 


(65) 


( 66 ) 


(67) 


Hence, after the transformation by J-expansion and truncat¬ 
ing the higher order corrections of order 0{t ~‘>^+^), we find 


K 


:= 


n[ 

n—1 


1 + g. 


-r 


d log t. 


f<N = fe¬ 


rn 


Our estimation protocol is almost the same as that for j3c- 
In the place of /<w, f>N is substituted because there is 
the matching region where the scaling behavior may be ob¬ 
served in />Ar, and then fc and unknown critical exponent 
qi will be estimated by utilizing extended PMS and self- 
consistent conditions (to be clarified later). 

From /;3> = ^ 

/^> = 1 - 12a; -f 104a;^ - lOOSa;^ -f 10416a;'^-, (69) 


and the (5-expansion at the expansion order N transforms the 
above result to fp^N := Disi[fp>N], 

fp>N = 1-12^ ^ ^ 

where the last term should be of the order . We note that, 
while the highest order of m is 25th, the highest order of 
fp>N is 24th, because of the cancellation of t in the numer¬ 
ator and the denominator in . 

The susceptibility at high temperature is also written with 
1/M = X (<C 1). First remind that %> = 1 -f 6^ -I- 30/3^ + 
0{I3^) 12^ . Then, by the substitution of {x) given in (l24l l 
into /3, we obtain 


f-fl04 


f^-1008 


f3+. 


',(70) 


X> = l + 6x-6x^ + 36a;3 - 270a;^ -f 2268a;^ 

-20436a;® -f 193176a;^ - 1890462a;® -f 18990892a;® 
-194709708a;^° -f 2029271688a;^^ 
-21435300372a;^® -f 228983179752a;^® 
-2469626018184x^4 -f 26855777435248x4® 
-294145354348974x4® -f 3242105906258220x44' 
-35935261094616124x4® 
-f400295059578038760x4® 


-4479014443566807276x®° 

-f50319506857313420376x®4 

-567383767790459777016X®® 

-f6418899321986117552400x®® 

-72838651914163555355012x®4 

-f828839976149614386374184x®®-, (71) 

and from/x> = x>'’/x>^ 

= 6x - 48x® -f 432x® - 4167x4 -f 42336x®-.(72) 


We thus arrive at 


/x>^=6 , hr-48 ^ x^+432 


• • •. (73) 


where 79 tv [/x>tv] := fx>N- For the sake of notational sim¬ 
plicity, let us omit the subscript N in //3(x)>(<)tv hence¬ 
forth. 


4.2 Biased estimation of v 

First of all, we like to confirm that our approach provides 
results consistent with (l23t . ( l59b and (l60l l by the use of the 
value uj = 0.84(4) in (ISFt . giving 

(?! = I = 0.42(2). (74) 

Such a biased estimation strengthen the validity of the ap¬ 
proach when successful and suggests possible improvement 
of the protocol in an unbiased estimation which will be ex¬ 
plored later. 

Before the estimation, we survey the behaviors of fpy, 
/^> and their derivatives. We find from Fig. 11 (a), (b) that 
the functions to the third order derivatives may exhibit scal¬ 
ings. The plots of log((—1)'“/^^^) and log((—1)^+4jW) 
shown in Fig. 12(a), (b) provide us more detailed informa¬ 
tion on the scaling levels. Both for /^> and f^y, the deriva¬ 
tives to the second order show scalings while the third order 
derivative is not sufficient. As for the fourth order one, we 
observe no linear-like behavior. This implies that the deriva¬ 
tives to the second or at most third order are appropriate for 
the estimation. The level of the scaling shown in the deriva¬ 
tives is less-qualified to that of the square model (see Fig. 
3(a)). 

As the first protocol (i), we consider the estimation of ix 
via the function fp^ only. When the first order ansatz fp^ = 

— 1/(2^) -f const X x~'^^ is used, the FDF to be considered 
is = [1 -f qi^{d/dlogt)]fp^ = —l/{2ix). Then, sub¬ 
stituting qi = ijj/2 = 0.84(4)/2, we search the stationary 
point around t = 0.11 ^ 0.13 where the scalings of fp^ 
and are observed (see the plot of (fc = 0,1) in Fig. 

11(a).) Then the stationary point at each order of = 21 to 
24 gives = 0.6078(34), 0.6077(35), 0.6095(33), 0.6094(33), 
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Fig. 11 (Color online) The graphs (a) and (b) show plots of and 
/^> at respective highest orders (N = 24 for and 25 for fx>) 
and their derivatives with respect to d/d log f to the 4th order. The 
gray lines indicate —l/{2v) for v = 0.6301 in (a) and 7/(2t^) for 
7 = 1.2373 and v = 0.6301 in (b). 



Fig. 12 (Color online) The graphs (a) and (b) show plots of 
log((—1)'“/^^^) and log((—1)^+^/^^) at respective highest orders 
{N = 24 for /j 9 > and 25 for /x>) for fc = 1, 2, 3, 4. 


respectively but these are not good approximation. Satisfac¬ 
tory estimation comes from the second order LDE, where 
the second order ansatz reads = —l/(2t^) + const x 
x~'^^ + const X a;“® . The second order LDE is given by 


'^I3< = 


1 + 9i 


-1 


dlogfJ - 


1 + 92 


-1 


dlogf 


/,< = -^.(75) 


Since l/(2i/) in the right-hand-side is nothing but q 2 , we can 
write the above LDE as 


^/3< = //3< + ( 9 l ^ + 92 ^) 4 < + (9192) V;3< = - 92 -( 76 ) 


By this installation of the self-consistent point of view, we 
can reduce the number of equations and the order of deriva¬ 
tives included in the estimation. We thus estimate v by the 
both use of (l76l l and the PMS condition, 

4< = f7< + + 72'')4< + (7i92)“7g = 0. (77) 

In the successive steps, we replace by . It is cum¬ 
bersome to remark this in each and every cases, we like to 
write LDEs only with and . We first solve the LDE 

= 0, and obtain involved in the correction as the 
function of t, the estimation pout. The solution denoted by 
<32/3 is given by 


<32/3 (t) ^ = 


/g + (2/u;)/g 
/g + (2/u;)/g' 


(78) 


Then, substituting the solution Q 213 (t) into <72 in both sides 
of (|26l), we look for the self-consistent point around, say, 
t = 0.11 ~ 0.12. The result is summarized in Table 4 in 
the ni row. Actually, at odd orders, the solution t of (l76l) just 
appeared from 23rd order. Eor lower odd orders, two curves 
of and —Q 2 p{t) have no common point and we have 
picked up approximate solution of t where ^/3> — (—<32/3 (f)) 
becomes locally minimum in the absolute value. Also at 
even orders, we consider that reliable solution just emerge 
from 22nd order, since from that order the second order 
derivative begins to show the tendency toward 

zero, after the peak at lower t, around the self-consistent 
point. We see from Table 4 that the result manifests good 
accuracy at high orders. We however point out one draw¬ 
back contained in the approach. That is, the scaling of Q 2/3 
which is deduced from the differentiation of (IbTT i that 

Q2/3(f)-' = - + <3(f-®+«^), (79) 

92 


is not seen even at higher orders such as (V = 23,24. As 
seen in Eig. 13 (a), (b), the curves of —Q 2 P show no sign of 
approaching to — l/{2v) represented by the dashed line. 

Now, to explore the possibility of handling series for the 
estimation, let us consider the use of to replace Q 2 P by 
a better function. The series /^< is expected to have almost 
the same structure with //3<; Though differences remain in 
the limit x ^ 00 and the amplitudes of each corrections, the 



























































14 


Hirofumi Yamada 




Fig. 13 (Color online) Plots of ^/ 3 > | (t) (blue) at = 23 and 24 

and Q 2 p{t) (red) for uj = 0.84. Horizontal dashed lines indicate the 
standard value of—l/(2iz) = —1/(2 x 0.6301). 



Fig. 14 (Color online) The graphs show plots of Q2x(^) (blue) and 
Q 2 ^(f) (Green) for u = 0.84 at respective highest orders {N = 25 
for the first and 24 for the second functions). 


both functions have corrections etc. The 

function /x> has an advantage, compared to that the 
derivatives show better behaviors than those of /^> deriva¬ 
tives (see Fig. 11 and Fig. 12). For instance, from Fig. 12, 
we find that (fc = 1,2) exhibits a bit longer linear-like 
behavior than that of . Thus, from the stationarity con¬ 
dition for 

= fx> + + <l2^)fx> + (9192)"^/x> = 0, (80) 


we obtain 

Q2x{t)~^ 




X> 


7(2) 

Jx> 


+ (2/tv’)/x> 

+ (2/w)/^f 


(81) 


which represents at the estimation point t. To demon¬ 
strate the differences between Q 2 i 3 {t) and Q 2 x{t) clearer, 
we plotted them at the respective highest orders in Fig. 14. 
The function Q 2 x{t) exhibits the rough scaling and the sta¬ 
tionary value at the top of the hill signals l/(2rz)-value (As 
Q 2 P obeys the scaling (l79l l, Q 2 X also obeys the same scal¬ 
ing. If the scaling could be of high grade, it should develop 
a plateau and the stationary point indicates q^^)- Our recipe, 
which is the point in this protocol (ii), is to replace (t) 
by Q 2 ^{t) in ^/3>. Here we note that for the combined use 
of and must be paid that the solution of 

= 0 at (iV -f l)th order should be substituted into A^th 
order ^/3>, since both functions are derived from {N + l)th 
order expansion of I3{M). By the substitution of Q2x(t) into 
(?2 in (f76l) . we search the solution obeying self-consistency 
in exact or approximate manner. For instance, as shown in 
Fig. 15 (a), the exact solution exists at At = 23 as the in¬ 
tersection of two curves. On the other hand. At = 24 case 
satisfies an approximate one, — {—Q2xit)) ~ 0(7^ 0) 
(see the behaviors around t ~ 0.12 plotted in Fig. 15(b)). 
The result is recorded in Table 4 as a r'2-sequence. As in the 
first protocol, the result is good at higher orders. We add one 
point to be noted that in this protocol the stationarity of 
at the estimation point is not exactly realized due to the re¬ 
placement of by Q 2 ^ . However, the local gradient at 
the estimation point is very small and the stationarity is still 
respected precisely. 

As the last protocol (iii), we consider the sole use of the 
PMS condition in the combination of and fp. The recipe 
uses the two conditions, 3^^ = 0 and ^ 0. The first 
condition gives Q 2 X (see (ISTT) ) and then its substitution into 
yields the solution t = t pms where the first derivative 
becomes zero or nearly zero. Then we obtain the estimate of 
V from (1751) . 

^/3>|iPMS =(82) 

The proper solution appeared just at 23rd order (non-proper 
solution at the order, giving = 0.62787, has omitted in 
Table 4). At this order, we emphasize that the self-consistency 
is also satisfied approximately, since 1/(2(32 x(^pms)) = 
0.62924 which is close to v = 0.63020 estimated by dS^ . 
This approximate realization of the self-consistent feature is 
observed for the first time at 23rd order. Though proper so¬ 
lution is not obtained at 24th order, it would emerge at least 
a few larger even order. 

We have plotted Vi {i = 1,2, 3) against the order N in 
Fig. 16. Though the extended PMS protocol (iii) needs a few 
higher orders for the assurance of convergence issue, we can 
say that all three protocols suggest the trend of the conver¬ 
gence to the correct limit. For the protocol (iii), however, 
the 25th order 1/M expansion (24th order for fp) is a bit 
too short. To summarize our estimation, we therefore focus 
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Table 4 oj-biased estimates of v via three protocols (i), (ii) and (iii): (i) estimates via = 0 and = — 92 - (ii) estimates U 2 via 

= 0 and = —<12 - (iii) estimates 1/3 via the PMS conditions = 0 and = 0. Parenthesis denotes the variation range of 
estimated v under the change of lo from 0.80 to 0 . 88 . 


order 

21 

22 

23 

24 

IZl 

0.62685(54) 

0.63004(63) 

0.63009(62) 

0.63013(62) 


0.62977(62) 

0.62929(54) 

0.63023(59) 

0.63007(47) 

JZ3 

0.62791(67) 

0.62787(76) 

0.63020(45) 

0.62782(78) 



1,0 III I. I. I. I i , , I . 

0.00 0.05 0.10 0.15 0.20 0.25 


Fig. 15 (Color online) Plots of ^/ 3 >|Q 2 ^(t) (blue) at AT = 23 and 
24 and Q 2 x{t) (red) for u! = 0.84. Horizontal dashed lines indicate 
the standard value of —l/{2v) = —1/(2 x 0.6301). In the last plot, 
we see that ^/ 3 >|Q,^ 2 (t) eome close with each others 

around t ~ 0.118 and approximately satisfies self-consistency repre¬ 
sented in <761 . 


on the self-consistent two protocols. Rounding the two es¬ 
timates of tz at the respective highest orders in protocols (i) 
and (ii) off to four decimal places, we hnd the both results 
agree with each other and yield single estimate. 


tz = 0.6301(6). (83) 

Here the number in the parenthesis implies approximate de¬ 
viation when uj = 0.80 and 0.88 are used. Note that the in¬ 
dicated range comes from the uncertainty of ui and not from 
some statistical origin. To compare our result (l83t with the 
reference value i/ = 0.6301(4) quoted in Ref. ITOl . our esti¬ 
mate is in excellent agreement. 



Fig. 16 (Color online) Plot of cj-biased estimates of iz via three proto¬ 
cols (i), (ii) and (iii): (i) Red plot with crosses indicates the result ui 
via = 0 and = —92 . (ii) Blue plot with boxes indicates 
the result U 2 via = 0 and = —92 ■ (iii) Black plot with 
triangles indicates the result 1/3 via the PMS conditions 3^'^ =0 and 
= 0 . Error bar denotes the variation range of estimated u under 
the change of ui from 0.80 to 0.88 (Lower and upper ends represent 
the limit for uj = 0.80 and 0.88, respectively). Dashed line indicates 
u = 0.6301 and the covered light gray zone indicates the range sug¬ 
gested in ED. 

4.3 Biased estimation of rj and 7 

The exponent p is estimated by using the first order ansatz 
= 'yl{2i/) + const x and the second order one, 
fx< = 'yl{2i') -f const x x~‘^^ + const x x~'^^. The re¬ 
sult from the hrst order LDE manipulated in the same way 
as //3> is not acceptable, however, since the estimate gives 
negative value of rj in some cases. Therefore, we need sec¬ 
ond order ansatz satisfying 

^x> = fx> + (9i ^ + '?2 ^)/x> + ('liQ'a) ^fx> = 

There are essentially two prescriptions to estimate p. 
One is, to avoid the bias as long as possible, to treat 92 as an 
adjustable parameter for the PMS. The PMS conditions to be 
used reads = 0 and ~ 0. From the first condition, 
we have Q2x(.t)~^ in exact agreement with (ISTI) and replace 

— / o') 

92 by Q2x{t) IB the second derivative 3^x>- search 

its zero or closest point to zero in the region f = 0.11 ^ 0.12 
where the emergence of the scaling is expected (See Fig. 
11(b)). By substituting the solution Q 2 x{t) into we ob¬ 
tain the estimate of 7/(2zz) and then 91, the estimate of 9. 
The other recipe is just to use estimated 92(= l/2zz) oh- 
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Table 5 oj-biased estimation of rj with the second order LDE. tji means the result from extended PMS and ri 2 from the simple substitution of 1^2 
and PMS. See the main text in detail. 


order 

22 

23 

24 

25 

771 

0.03367(61) 

0.03586(44) 

0.03672(31) 

0.03695(25) 

V2 

0.03770(50) 

0.03755(48) 

0.03761(47) 

0.03758(46) 



Fig. 17 (Color online) Plots of tu-biased estimates of rj via two pro¬ 
tocols: Black plot indicates the results via PMS conditions = 0 
and ^ ~ 0 without using estimated u. Blue plot indicates the re¬ 
sults via ^ = 0 and with the substitution of U 2 . EiTor bar denotes 
the variation range of estimated rj under the change of u! from 0.80 
to 0.88 (Lower and upper ends represent the limit for uj = 0.80 and 
0.88, respectively). Dashed line indicates rj = 0.0364 and the covered 
light gray zone indicates the range suggested in 1591 . 

tained in the previous subsection into and search the 
least variation point of at the appropriate region. The 
point provides the estimate of 7/(2tz) and 772 of 77-estimate. 
As an representative value of <72 in this prescription, we used 
1^2 in Table 4 . The result of the two estimations is summa¬ 
rized in Table 5 and Fig. 17 . At the highest two orders, all 
recipes give a bit larger values compared with the standard 
one, (I59I 1. 

As the last biased estimation task, we investigate the 7 
estimation. First we like to comment that the combination 
fxlf /3 '■= f-y gives —7 in the M —>■ 0 limit and, if 
would allow us effective ^-expansion, it would be a best 
function for our purpose. However, has a pole at a; = 
0.2388637 • • •, which is confirmed by its diagonal Fade ap- 
proximants to high accuracy and this spoils the power of 
the ( 5 -expansion. Thus, we use the estimated jz and 77 to ob¬ 
tain 7 through 7 = 2v{l — 77/2). In this protocol, two op¬ 
tions are possible: One is to substitute 771 and the other is 
to substitute 772. Thus, we have two candidates of estimates, 
7 i = 27/2(1 - 771/2) and 72 = 27/2(1 - 772/2). The result 
of estimates are summarized in Table 6 and Fig. 18 . Though 
the accurate estimate is not achieved yet, the trend of con¬ 
vergence is observed and the results at higher order implies 
agreement with the reference value dhOl l. 

To summarize this subsection, we hnd rough agreement 
with the reference results ( | 59 ] | and dhOl) . In the estimation of 



Fig. 18 (Color online) Plots of 7<j-biased estimates of 7 via two pro¬ 
tocols: Black plot indicates the results via PMS conditions = 0 
and ' ~ 0 without using estimated !/. Blue plot indicates the result 
under the simple substitution (72 = l/(2i/2) with PMS condition on 
^x>- protocols, we used Fisher relation 7 = 2i/(l — 77/2) to 

compute 7 with 1/2 ■ Error bar denotes the variation range of estimated 
7 under the change of u from 0.80 to 0.88 (Lower and upper ends 
represents the limit for cj = 0.80 and 0.88, respectively). Dashed line 
indicates 7 = 1.2372 and the covered light gray zone indicates the 
range suggested in 1601 . 

both exponents, it seems that the PMS protocol is slightly 
superior to the simple substitution of pre-estimated 7 /. When 
(72 does not obey in the included function some constraint, 
such as the self consistent one, it maybe more suitable to be 
dealt with as an adjustable parameter,. 

4.4 Unbiased estimation 

Here we attempt to estimate critical exponents in unbiased 
and self-contained manner. We first address to the estimation 
of uj and 7 /. As many terms of the corrections are included 
in the estimation becomes accurate. However, to fix 

associated unknown exponents, derivatives of target func¬ 
tions should be involved as many and then large order series 
in 1/M expansion is required. Confined with 25th order, it 
is safe to rely upon the two parameter ansatze, 

// 3 < = — ^—h const ■ + const ■ , (85) 

/y< = — -f const ■ x~'^^ + const ■ x~'^^. (86) 

^ 2 v 

As the hrst attempt, let us discuss the extended PMS us¬ 
ing conditions = 0 which gives ( 3 iy(f) and 

Q 2 x{t) for the corresponding exponents <71 and ( 72 . Then, 
substituting them into the last condition that becomes 
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Table 6 Lj-biased estimation of 7 with the second order LDE. In evaluation of 7 at order N through j = 2u{l — ril2), we used u = U 2 estimated 
in the previous analysis at order + 1. 71 = 2i/2(l — tJi/S) and 72 = 2i/2(l — t?2/2). 


order 

22 

23 

24 

25 

71 

1.23833(108) 

1.23602(94) 

1.23732(87) 

1.23686(100) 

72 

1.23579(91) 

1.23496(76) 

1.23676(86) 

1.23646(64) 


least sensitive at < = tpMS in the appropriate region, we 
obtain estimates of qi by Qix{ipMs) = Q* (* = 1)2). Since 
qi =a;/ 2 andg 2 = l/( 2 t 2 ), may give estimates 

of oj and i/. However, from an experience in Pc estimation, 

we have understood that i/pMS := ^l3>\t=tpMs,gi=QixitPMs) 

is more reliable for ly estimate (see (fTSll ). The result reads, 
for example at At = 24, vpMS = 0.63537, v*{= l/SqJ) = 
0.41331 and u}*{= 2ql) = 1.37326. The results at other 
orders are also of the same level or worse. The problem is 
found in the large discrepancy between the two estimates of 
V by vpMS and v* (former is better than the later) and this 
would have some connection to largeness of estimated uj. 
The circumstance is quite similar when = 0 is 

used instead of = 0. The pure PMS protocol is 

therefore not acceptable. 

As explicitly confirmed in the previous biased compu¬ 
tation, one of the satisfactory features to be possessed in 
good protocols is an approximate realization of the self- 
consistency. From this point of view, we consider that the 
best way is to employ the second protocol (ii) presented in 
the former subsection, which protocol is a hybrid of using 
PMS and the self-consistency in the combination of and 

fp- 

Let us summarize here the conditions of the protocol (ii)- 
unbiased as follows; 

- The self-consistency for the parameter <72 is rigorously 
respected at the estimation point. 

- Use the stationarity condition = 0 which gives 
92 ^ Q 2 x{t) ^s its effective function of estimation point 
t. 

- The upper limit of the observed scaling behavior of Q 2 X (0 
is matched with the self-consistent point. 

According to the above protocol, the estimation goes as 
follows: For example, consider the case at A* = 24, the high¬ 
est order of (pp. Taking three conditions into account, we 
impose the coupled equations 

^P>\q 2 =Q 2 x = “Qzx, (87) 

gg = 0. (88) 

The second equation comes from the interpretation that the 
upper limit of the scaling of g2x would be attained at the 
stationary point. We then obtain the solutions, 
t = 0.11793799038 ■ ■ ■ = tsc andoj = 0.7451128251 • • • = 
oJsc- Then, the substitution of the solution into ^ 0 >\q 2 ^ 



l.Ol " I I' I' I —^^^^^^^^^^^^^^^^—7 

0.00 0.05 0.10 0.15 0.20 0.25 


Fig. 19 (Color online) and —Q 2 x{t) as functions of t 

which specifies the solution of stationailty condition = 0. Com¬ 
pare with plots in FIG. 13(c),(d). Horizontal dashed lines indicate the 
standard value of—l/(2iz) = —1/(2 x 0.6301). The common point 
of the two curves indicates the acceptable self-consistent solution. 


or g 2 x yields v = l/{2q2) = 0.628830753 •• • = vsc- 
Next, consider the case at A* = 23. At this order the set 
of solutions is absent and we loosen the second condition 
such that g2x should approximately satisfy the stationar¬ 
ity at the self-consistent point t = tsc\ That is, we in- 
stead use Q^x = 0 at f = tsc (Above 20th orders, we 
found that at even orders tsc is obtained as the station¬ 
ary point of Q 2 x and at odd orders as the least variation 
point in the scaling region.). This modification then gives 

tsc = 0.11673398933-and wsc = 0.7120047317 
Thus, we obtain vsc = 0.6278846723 ■ • -. For the graphi¬ 
cal representation of the situation, see Fig. 19 and Fig. 20. 

The results from 20th to 24th order are summarized in 
Table 7 and Fig. 21. We observe that both of ijJsc and vsc 
are monotonically increasing with the order N. For uj re¬ 
sults, this may explain the smallness of use such that the 
incorporation of un-computed higher order terms would lift 
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Fig. 20 (Color online) ~Q2xi^) if vicinity of 

the self-consistent solution represented by the cross mark at 24th order. 
Horizontal dashed line indicates the standard value of —l/(2u) = 
-1/(2 X 0.6301). 


Table 7 Unbiased estimation of lo and u via the self-consistent PMS 
protocol (\\)-unbiased. 


order 

20 

21 

22 

23 

24 

use 

Use 

0.31643 

0.60424 

0.41711 

0.61507 

0.61357 

0.62526 

0.71200 

0.62788 

0.74513 

0.62883 



20 21 22 23 24 

Fig. 21 (Color online) Plots of estimated u! (plot (a)) and i/ (plot (b)) 
via the protocol (ii)-unbiased. 


up the higher order estimates near (l58l l. The estimate of ly at 
the present highest order is slightly smaller than (l 2 Tt but the 
sequence {lysc} indicates the value ~ 0.63 and essentially 
consistent in the world average. 


Table 8 Estimation of p and 7 biased by use obtained in the self- 
consistent PMS protocol (ii)-unbiased. Lists of two kinds of 7 come 
from the following : Upper one denoted by 'ysc shows the result ob¬ 
tained by use of the estimated use and the lower one represented by 
the asterisk by use of u = 0.6301 quoted in d231 from ref. P). The 
later sequence is added to guide how estimate of 7 is affected by u- 
value. 


order 

21 

22 

23 

24 

25 

Vse 

0.01186 

0.02474 

0.03355 

0.03578 

0.03633 

7se 

1.20131 

1.21869 

1.22954 

1.23387 

1.23481 

7’ 

1.25272 

1.24461 

1.23906 

1.23766 

1.23731 


Next, we turn to estimate 77 and then 7 . We use use 
just obtained in the protocol {n)-unbiased. The estimate of 
V is not required since it is desirable to avoid as possible 
the bias from the previous estimation. Then, we use dSTT i to 
substitute it into the second derivative of and select the 
point satisfying 

~ 0 ( 89 ) 

where ^ means the exact zero (at odd orders) or approxi¬ 
mate zero (even orders). Having obtained the solution fpMS, 
we simply substitute it into and identify 


2v 2 


q2—Q2x{^PMs)yt—ipMS 


(90) 


We obtained the results shown in Table 8 and Fig. 22 (a). 
The estimate of 7 is simply done through 7 = 2u{1-t]/2). 
See the result in Table 8 and Fig. 22 (b). We first remark 
that, at iV = 25, tpMS exactly agrees with tsc obtained 
in the u estimation. This is because the second stationarity 
condition = 0 is equivalent with “ 9 (see (l 88 l l). 

When ^^4. 7 ^ 0 at tpMS, tpMS shows just a slight differ¬ 
ence from tsc- These facts mean that the value of Q 2 X at 
tpMS yields exactly the same or very close value to use at 
respective orders (Remind that Q 2 X converges to 1/ {2u) in 
the scaling limit). The protocol (ii)-unbiased thus fulhlls the 
unihed treatise of estimations of various exponents. 

The estimate of 77 at 25th order agrees excellently with 
the quoted value (l59l l IfTOl . The lists of the 7 estimation are 
composed by two sequences. The upper sequence stands for 
the results by use of use shown in Table 7. The lower se¬ 
quence asterisked shows the results by use of u quoted in 
(l23ll as the standard lITOl . The result of 7 at = 25 is 
slightly lower than the standard one dhOl l though the trend 
of the sequence implies acceptable limit. The list of 7 in the 
upper sequence is biased from use and it would be inter¬ 
esting to see that the use of (| 2 ^ gives excellent agreement 
with the standard 7 quoted in (I 6 OI 1 . This indicates that our 
approach to the critical exponent estimation is totally satis¬ 
factory. 
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0.05 

0.04 

0.03 
Vsc 
0.02 

0.01 

0.00 
1.28 

1.26 

1.24 

1.22 

1.20 

1.18 

Fig. 22 (Color online) Plots of estimated rj (plot (a)) and 7 (plot 
(b)) biased by ujsc obtained in the self-consistent PMS protocol (ii)- 
Linbiased applied io Red cross represents the data 75 c ^nd blue 
triangle the data 7*. 
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by the numerical plots of relevant functions. This feature en¬ 
ables us to identify which solution among the ones produced 
by the extended PMS or self-consistency condition should 
be really relied upon. 

The works of Guttmann, Butera and Comi ||5]|6l, using 
solely the high temperature expansion with the aid of dif¬ 
ferential approximants method, used shorter series and the 
direct comparison is not possible. We just mention that our 
estimates of and 7 are better than their results. In jS], 7 
is estimated as 1.2431(24) and in lb), v = 0.634(2) and 
7 = 1.244(3). However, their estimates of Pc is better than 
ours; Pc = 0.221657(7) in Q and 0.221663(9) in El. In 
fact, we have expected that the improved estimate of critical 
exponents may help to improve the pc estimation. However, 
we have not obtained yet a solid scheme for the essential 
improvement. 

Finally, we remind that the estimation protocol (i\)-unbiased 
fulfills the conditions of self-consistency in exactly and 
the approximate (or exact in cases) stationarity of Q 2 X and 

at the estimation point. These are also valid for func¬ 
tion for rj estimation. We note that those conditions should 
be respected when the functions under the investigation are 
really in the scaling region. The set of those conditions is 
strong and tightly bounds the estimation procedures. If one 
can find more flexible variant, the accuracy of estimates might 
be improved. 


5 Concluding remarks 

To summarize the paper, we have investigated the estimation 
of critical quantities in expansion in a mass argument. First, 
the critical inverse temperature Pc has estimated through 
LDE with constant coefficients related to the critical expo¬ 
nents. With the aid of the ^-expansion and the extended PMS 
condition, we have obtained at the highest order 25th, Pc = 
0.221687, larger than the standard upper limit ^ 0.22166. 
As a by product, we have also obtained the rough estimate 
of V at third order LDE. 

Independently, the improved estimations of v and other 
critical exponents have been attempted by employing P^'^'>/P^^'> 
and /x as the target functions. Under the hybrid pro¬ 
tocol of PMS and self-consistency point of view, we have 
attempted biased and un-biased estimations of critical expo¬ 
nents. Biased hy oj = 0.84(4), we obtained results in good 
agreement of world averages. In unbiased approach, satis¬ 
factory results are obtained for ly, rj and 7 in self-contained 
way; The best estimate of uj is 0.74513 and smaller than the 
range indicated in ref. m, though the obtained sequence 
may probably grow to the value around ^ 0.8 at 
large enough orders. Exponent v has estimated to be 0.62883, 
slightly smaller compared to (l23T l. In contrast, the estimate 
77 = 0.03633 agrees very well with ( l59l l. We emphasize that 
in our approach, the scaling behavior is directly observable 
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